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ABSTRACT 


The study of chlorophyll content in urban vegetation 
leaves provides important information about the 
plants health and ability for CO2 absorption. Recent 
advances in optical remote sensing led to the 
development of improved methodologies for remote 
estimation of chlorophyll content. For instance, the 
hyperspectral index NAOC (Normalized Area Over 
reflectance Curve), proved to be a promising 
indicator. 

In this work we studied tree chlorophyll content in 
the city of Valencia (Spain) based on measurement 
of chlorophyll taken from four representative tree 
species: plane tree (Platanus x. acerifolia), canarian 
palm (Phoenix canariensis), hackberry  (Celtis 
australis) and mulberry (Morus alba). Measurements 
were made during summer 2011 within the 
European project BIOHYPE (BlIOmonitoring of 
urban habitat quality with airborne HYPErspectral 
observation). For 320 leaf samples, chlorophyll 
content were measured both in the laboratory and by 
using a SPAD-502 chlorophyll meter. These data 
were correlated leading to a best fit with power type 
functions, with coefficients of determination R2 
above 0.86. During the field campaign an aircraft 
flew over the city with onboard a hyperspectral 
sensor CASI (Compact Airborne Spectral Imager) 
and obtained several images with a resolution of 1 m. 
From the CASI data we first calculated the NAOC 
index and secondly related 1t to the laboratory 
chlorophyll content measures. However, in order to 
adapt the index to remotely estimated chlorophyll 
content over all urban vegetation, extra data of low 
growing vegetation coming from earlier campaigns 
was added. It led to a high resolution chlorophyll 
content map over the study area of Valencia. 


Keywords: Urban trees, chlorophyll, NAOC, SPAD 


1. INTRODUCTION 


Chlorophyll is a natural pigment found in green 
plants, and its content 1s an important variable in 
assessing the plants physiological state and 
photosynthetic performance, the basic process for the 
vegetation life(Wang and Li, 2012; le Maire et al., 
2008). Vegetation in urban areas, particularly trees, 
provides immense benefits to urban environment 
such as the main absorbers of CO». The 


quantification of Chl from urban vegetation allows us 
not only to gain insight in the plants” health status 
but also in their ability to regenerate the urban 
environment (McPherson et al., 2011). 

However, traditionally the measurement of Chl is a 
laborious and destructive process that requires 
sampling and chlorophyll extraction in a particular 
area Of the leaf by dissolving chemicals in the 
laboratory. Moreover, Chl concentration 1s known to 
be heterogeneously distributed across different leaves 
of the same plant, depending on sunlight exposure 
and the size and age of the leaves. And even 
within a single leaf the Chl concentration is not 
uniformly distributed. Therefore, in order to infer Chl 
in a particular tree or a large area of vegetation 1t 
requires the collection of large number of 
measurements, which cannot be performed in the 
laboratory for practical reasons. Alternatively, to 
overcome this limitation, instruments have been 
designed that are capable non-destructive and fast 
estimating of Chl, directly by measuring the leaf light 
absorption at certain wavelengths. These instruments, 
called chlorophyll meters, provide a numerical 
value which is related to the actual Chl through a 
calibration process that may depend on the 
vegetation type and its phenology. One of the widely 
used chlorophyll meters is the SPAD-502 and is used 
1n this study. 


However the study of Chl over large areas such as 
over a city can be quite challenging. Here, an 
attractive alternative 1s to make use of optical 
remotely sensed images (Gamon «¿€ Qiu, 1999). 
During the last decades various techniques have been 
developed that estimate Chl from images acquired 
from sensors mounted on an air- or spacecrafts 
(Wang á£ Li, 2012). These techniques are usually 
based on the characteristic chlorophyll spectrum with 
two peaks in the blue and red visible part. However, 
also other plants pigments absorb light in the blue 
region. This 1s why other optical Chl estimation 
techniques prefer to make use of bands located 
between 600 and 700 nm, which is the spectral region 
where chlorophyll is the main absorber. Thanks to the 
development of hyperspectral sensors, which can 
obtain a continuous spectrum in this region, or new 
superspectral sensors such as ESA's forthcoming 
Sentinel-2 mission, which also incorporates multiple 
bands in this region (Delegido et al., 2011b), new and 
more powerful Chl detection methods are being 
developed. One such example is the NAOC 
(Normalized Area Over reflectance Curve) index, a 
spectral index that was initially developed for Chl 
mapping over heterogeneous crops areas (Delegido et 
al., 2010). 

This brings us to the aim of this study, which is to 
recalibrate the NAOC index so that it enables Chl 
estimation over urban trees. To do so, four most 
representative trees species in the Valencia city 
(Spain) were first sampled by using SPAD-502 
measurements and through destructive Chl extraction 
used for calibration. The Chl  estimates were 
subsequently linked to NAOC and then applied to an 
alrborne image for upscaling Chl mapping over 
Valencia city. 


2. MATERIAL AND METHODS 


2.1. Test site 

The data used in this study were obtained within the 
BIOHYPE project (BlOmonitoring of urban habitat 
quality with airborne HYPErspectral observation) 
which has the objective of developing, testing and 
validating a passive biomonitoring methodology 
based on hyperspectral observations from an aircraft 
level and with an experimental field approach. 

Within the BIOHYPE project there have been 
different field campaigns organized in order to 


obtain a series of measures of different biophysical 
parameters, in addition to their radiometric and 
hyperspectral imaging obtained by a sensor airborne. 
One of these campaigns took place in Valencia from 
May to October 2011. This work was carried out 
with data obtained during the Valencia campaign. 
Valencia (39%28' N, 09222" W, 15 m a.s.1.) is the third 
largest city by population in Spain, with 798,033 
inhabitants (INE, 2012), and an area of 134.65 km”. 
It 1s located on the Mediterranean coast of the Iberian 
Peninsula, on the banks of the Turia River. The 
climate is typical Mediterranean with mild, wet 
winters and hot dry summers, and a yearly-average 
temperature around 18 *C. Seasonal rainfall is 
minimal in summer and maximal in autumn and 
spring. The city is divided by a large park called 
Gardens of Turia which are built in the old Turia 
riverbed. Another large park, Viveros, 1s closely 
situated near the Gardens of Turia (Figure 1). For the 
whole city of Valencia 116,745 trees are counted in 
May 2010. Tree selection was made based on the 
species distribution maps provided by the City 
Garden Council Service, taking into account the 
distribution and occurrence of different species and 
leaf characteristics. Four dominant tree species were 
selected: London plane (Platanus x. acerifolia 
(Aiton) Willd), the Canary Island date palm (Phoenix 
canariensis Chabaud), European nettle tree (Celtis 
australis L.), and White mulberry (Morus alba 
L.). The city counts respectively 8.922, 3.513, 
2.935 and 2.193 trees of P. x acerifolia, M. 
alba, C. australis and P. canariensis, presenting in 
total more than15% of all city trees. 

In order to sample trees growing, ten city locations 
were chosen with different environmental conditions. 
In total 40 trees were sampled, ten for every species. 
Twelve trees at three of those ten locations were 
selected in order to take full crown measurements by 
crane at the low, middle and top layer of the crown. 
These three sites (Turia, Viveros and Glorieta) are 
three parks located relatively close to each other 
(Figure 1), but differ in area and traffic intensity 
surrounding them. The trees at Glorieta are 
surrounded by high buildings, Viveros, on the other 
hand, is one of the largest city parks where trees are 
somehow more isolated from traffic, and Turia park 
1s the old riverbed of the river Turia and now 
planted with trees, forming an important recreation 
zone of the city. 





Figure 1. Test site location. 1. Península Ibérica; 2. Valencia (with flight line); 
3. a. Viveros; b. Turia; c. Glorieta 


2.2 Chlorophyll extraction 

From each of the twelve trees a branch was cut at 
each of the three canopy levels (bottom, middle, top). 
These branches were immediately transferred to the 
laboratory in buckets filled with water. In the 
laboratory 3 leaves of each branch were chosen. 
Relative chlorophyll measurements were made with 
SPAD-502 (Konica Minolta Optics Inc., Japan) 
before samples were cut with a metal cylinder and 
stored in a freezer at -80% C. The samples were cut 
into the center of the distance between the stem and 
the apex of the leaf, and between the edge of the 
blade and median nerve. Upon chlorophyll extraction, 
samples were drawn from the freezer to measure the 
chlorophyll content by spectrophotometry in the 
laboratory. Each sample was ground in a porcelain 
mortar with 5 ml acetone and kept cold in an ice bin. 
After the sample was centrifuged for 5 minutes at 
10000 rpm and then measurements were made with 
a spectrophotometer. The absorbances were measured 
according to Lichtenthaler and Buschmann (2001). 
Later, after an initial inspection of results 1t appeared 
that SPAD values were generally above 30 or 40 
ne*cm?, new measures were taking from 140 
leaves (35 of each tree species) to complete the solid 
curve for SPAD-Chl calibration. In total 320 samples 
were measured and used for calibration. 


2.3. SPAD calibration 

The SPAD instrument provides a relative value for 
chlorophyll content. The objective was to calibrate 
the SPAD with measurements from chlorophyll 
extraction in the laboratory. Several authors have 
addressed this problem and a wide variety of 
calibrations of SPAD for different plant species can 
be found in literature. Most of them provide 
calibrated linear, exponential or  polynomial 
equations.  Often, a linear  equation was 
encountered especiall y when chlorophyll 
concentration ranges are narrow (Richardson et al., 
2002). However, when Chl ranges are broader, the 
trend can be better described by an exponential or 
polynomial relationship. Overall, Chl values for tree 
species can vary between 0 and 150 ug*cm” 
(Marenco et al., 2009), while Chl values for crop 
species vary between O and 70 (Uddling et al., 2007; 
Delegido et al., 201 1a). 


2.4, Remote sensing data 

On 7 September 2011, a flight campaign was 
made by INTA (National Institute of Aerospace 
Technology) with the CASI sensors on board. 15001 
CASI (Compact Airborne Spectral Imager) is a 
hyperspectral sensor in the range of 350-1200 nm 
that can work in several acquisition modes. 


CASI images were obtained during 4 flightlines 
over different parts of Valencia covering the 10 
sample sites. For this work we used the CASI 
flightline that covered the three sites where 
measurements were conducted at three crown levels 
in the trees, configured with 144 spectral bands with 
2.4 nm and a pixel size of 1.6 m along the flight track 
and 1.0 m across. This is the best configuration for 
the NAOC calculation because providing a high 
bands number in the region of interest that allows 
calculating the reflectance integral and also provides 
high spatial resolution which enables distinguish 
individual trees. The image was atmospherically 
corrected by the method of Guanter et al. (2007) and 
georeferenced. 


During “image acquisition, two ground teams 
conducted spectral measurements and measured 
various biophysical parameters. A team with a mobile 
crane sampled Glorieta and Viveros, while another 
was stationed in Turia. 


2.5. Chl remote sensing estimation from NAOC 
The aim of this study was to extend the 
application range of the  hyperspectral index 
NAOC (Normalized Area Over reflectance Curve) 
for trees with high chlorophyll content. This index 
was initially developed for remotely detecting Chl 
over crop types and here 1t is evaluated whether this 
index enables remote estimation chlorophyll content 
over urban trees. 


The NAOC index is calculated by the formula (Delegido et al., 2010): 


b 
ao 


NAOC=1 - 


Pmax (b=a) 


(1) 


Where p is the reflectance, A the wavelength, Pmax 18 the maximum far-red reflectance, corresponding to reflectance 


at the wavelength “b”, and “a” and “b” are the integration limits surrounding the Chl well centered at -670 nm. In a 
previous work (Delegido et al., 2010), we obtained that this index reaches a maximum linear correlation with the 
chlorophyll content of leaves in a wide variety of crops, when the selected bands (a and b) are 643 and 795 nm. The 
method has been successfully validated against Chl data coming from a wide variety of crop types (Delegido et al., 
2011a). The NAOC index tends to be more sensitive to changes in Chl than other remote sensing indices, because 
of the relationship is more linear than other indices that have the saturation problem at high values. 


3. RESULTS AND DISCUSSION 

3.1. SPAD calibration 

The 320 chlorophyll values measured in the 
laboratory according to the method Lichtenthaler 
and Buschmann (2001) were plotted against the 


SPAD-502 recordings for all measurements in the 
four tree species in Figure 2. It led to a maximum 
correlation when using an exponential function 
according to: 
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Figure 2. Total chlorophyll measured in the laboratory as function of SPAD 
measurements of the four species selected with the fitted curve using equation 2 


When looking closer to Figure 2, 1t can be observed 
that the curve is marked by high Chl values 
exclusively coming from Phoenix (P. canariensis) 
leaves. It should be noted that Phoenix species are 
characterized by a specific physiology with a 
pronounced hardness and thickness of the leaf width. 
Hence, because of this distinct leaf physiology, a 
more correct approach would be to seek for 
correlations for each tree species independently. It 
was further noted that in each of these graphs the 
lowest Chl values, e.g. up to 35 SPAD adjusts well to 
a linear relationship, which may explain why many 
authors have proposed a linear fit especially for 
narrower ranges of Chl. It can be observed that both 
Celtis and Morus span a similar Chl range, from 0 to 
130 ug*cm”, while Platanus measurements never 
exceeded 70 ug*cm”. Remarkably, Phoenix Chl 
distribution appeared to be clustered towards two 
cloud points, with a first cloud for Chl <15 ug* cm”, 
and then a second cloud spanning values of 80 
ne*cm” to higher up until 200 ug*cm?. The Palm 
species behaved differently than the other species 
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because of being an evergreen tree with steady 
chlorophyll content throughout the seasons. Here, low 
Chl values came rather from leaves that are 
approaching the end of its vegetative period and start 
to senescent. 


3.2. Chlorophyll content mapping over Valencia 

A following step was to calculate the NAOC 
index from airborne spectral data and relating 1t to 
Chl in-situ measurements for Chl mapping of over 
the city of Valencia. Prior to be able retrieving 
actual chlorophyll content the NAOC index must be 
calibrated. To do so, we used with the reflectance 
spectra of the atmospherically corrected CASI 
image, from which we calculated NAOC. These 
values were  correlated with corresponding 
chlorophyll content measured in-situ at the 12 trees 
during the flight campaign (Figure 3a). For the 
measurement of chlorophyll from each tree, we took 
the arithmetic mean of all chlorophyll measures in the 
different sheets and calculated the standard deviation 
of each mean. 
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Figure 3. Chlorophyll content according NAOC. The bars indicate the standard deviation of the Chl 
measures. (a) All species of trees (b) Without Phoenix 


It appeared that the three canopy class values of 
Phoenix species deviated from the general trend. 
This might be due to the fact that Phoenix leaves are 
very narrow and sparsely organized on top of the 
trunk, so the fractional cover 1s very low for the 
majority of the crown area as seen from above. From 
a pixel point of view, this results into a mixture of 
tree crown and background for the pixels covering 


Chl (ugtem 5) = -234.76 + 619.32 *NAOC 


Palm trees. In addition, Phoenix leaves are much 
thicker than other tree leaves so underestimates the 
reflectance and chlorophyll content is difficult for a 
single spectral index present similar behavior in palm 
and other trees. For these reasons we have 
eliminated the three mean values of Palm in further 
processing (Fig. 3 b). The Fig. 3b data shows a linear 
trend that can be fitted to the following equation: 


(R? = 0.857) (3) 


Equation 3 would estimate the chlorophyll content 
for the great majority of trees in the Valencia city. 
However, this fit did neither include low 
vegetation (e.g., grass, herbs, bushes...) nor non- 
green targets (e.g., roads, buildings), so 1t would 
inevitably wrongly estimate chlorophyll content on 
those targets, e.g. when mapping Chl over a city. In 
order to resolve this situation, we have incorporated a 
calibration from previous ESA “s campaigns (SPARC, 
CEFLES2 and Sen3Exp, see details in Delegido et 


al., 2011a) over agricultural targets that show similar 
chlorophyll levels that those. It was found possible to 
find a single function relating the Chl and the NAOC 
will both low plant and trees, however the equation 
was no longer linear. Figure 4a shows the values of 
NAOC for all targets: trees, grass and low vegetation 
(similar to crops) and a non-vegetated surface (with 
Chl = 0). These values can be fitted according to an 
exponential equation: 


7.169 Z 


Chl (ugtem 5) = 67.74*NAOC + 5071.3*FNAOC (R” =0.,854) (4) 


To end with, by applying Eq. 4 on the CASI image 
a cChlorophyll content map has been generated 
(Figure 4b). In the map, trees with high chlorophyll 
content such as Celtis (about 100 ug* cm? in purple) 
can be clearly distinguished from the lowest content 
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4. CONCLUSION 


The measurements with the SPAD-302 meter over 
different tree species show a strong correlation with 
Chl obtained in the laboratory. Equation 2 provides 
the calibrated empirical relationship as obtained by 
four representative species trees of in the Valencia 
city and many other Mediterranean climate cities. 
The best performing correlation was obtained by an 
exponential-type function, but for lower values also a 
linear calibration would be acceptable. 





trees (Morus and Platanus in green and blue, over 60 
ne* cm). Light green and yellow areas represent low 
vegetation and orange-red correspond to  non- 
vegetation surfaces. 


(b) 
Figure 4. (a) Chlorophyll content according NAOC from crops and low vegetation, trees and bare soil. 
(b) Chlorophyll content map (in ug*cm) of test site 


Due to distinct physiological characteristics between 
Phoenix canariensis and other species (Platanus 
hispanica, Celtis australis and Morus alba), there 
were however clear differences observed in Chl. 

The leaf chlorophyll content of the palm tree was 
mostly sampled between 100 and 200 ug*cm * while 
Chl did not exceed 140 ug* cm “in the other species. 
The pronounced thickness of the palm leaves leads 
that both the SPAD measurements as the reflectance 
NAOC underestimated Chl. This requires a more in- 
depth study on these kinds of thick leaves, such as 
analysis on both the  reflectance as well 
transmissivity. 


As a next step, the hyperspectral NAOC index was 
calculated from the CASI sensor and successfully 
validated against the in situ measured Chl. Again 
there appeared to a pronounced different behavior 
between the palm and other species due to its 
aforementioned characteristic physiology. Besides, in 
contrast to the other 3 tree species, the specific shape, 
size and distribution of palm leaves causes that the 
fractional cover of the palm tree is low, which 
implies a relatively greater influence of soil 
background contamination. For this reason, the palm 
data have been removed from the fit between Chl and 
NAOC data. Alternatively, by  incorporating 
averaged data from various agricultural sites 
originating from earlier campaigns, a robust 
empirical relationship between Chl and NAOC 
estimates was obtained. This relationship can serve 
Chl mapping over low vegetation types (e.g., grass, 
shrubs.), bare soil as well trees. Finally, the 
developed method allowed us generating a Chl map 
over the Valencia city. 
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ABSTRACT 


Environment, climate and human activities play an 
important role in determining vector distribution 
and vector-borne diseases epidemiology. A better 
understanding of the contribution of those 
ecological factors on the spatio-temporal dynamics 
of mosquitoes may lead to improvement of 
controlling measures. In the present study, it 1s 
thought that the conceptual approach (CA) of 
Tele-epidemiology developed by the French 
Spatial Agency and its partners could be applied 
to spatio-temporal mapping of entomological 
malaria risk in urban settings and entomological 
dengue risk. In Dakar, Senegal, a large amount of 
entomological data was collected (2007 - 2010). 
Appropriate satellite environmental and 
meteorological data were acquired in order to 
model the presence of water bodies, the presence 
of Anopheles gambiae s.l. larvae in those water 
bodies and the An. gambiae s.l. Human Biting 
Rates. In Tartane, Martinique, entomological field 
data consisted in the routine surveillance data 
collection. Very high spatial resolution data were 
used to map the environmental elements that may 
be associated with the presence and type of Aedes 
aegypti breeding sites in houses. In Dakar, the first 
mapping level predicted the potential An. gambiae 
s.l. larval habitats (depending upon the predicted 
presence of water, the environmental and 
meteorological conditions). The second level was a 
map of the An. gambiae s.l. HBR (driven by 
surfaces and potential productivity of larval 
habitats, urbanization levels and rainfall amount). 


In Martinique, the basic information has been 
prepared in order to further apply the Tele- 
epidemiology CA to the mapping of dengue 
vectors. 


Keywords: Remote Sensing, malaria, dengue, 
entomology, risk maps 


1. INTRODUCTION 


Malaria is caused by a Plasmodium  parasite 
transmitted among humans by the female 
mosquitoes/vectors of the Anopheles genus. 
Dengue is an infectious disease caused by one of 
the four serotypes (DEN-1 to DEN-4) of a dengue 
virus, transmitted to humans by the bites of infected 
female mosquitoes of the genus Aedes. The location 
of the vectors larval habitats and their dynamics are 
the primary determinants of the spatial and temporal 
distribution of malaria transmission. Transmission 
units are defined as the systems where the larval 
habitats are the sources for the  pathogen 
transmission [Carter et al. (2000)]. Indeed, 
mosquitoes emerging from larval habitats have a 
somewhat limited flying range which depends on 
environmental conditions. Malaria and dengue risk 
1s thus heterogeneous in space and time, as it 1s 
driven by the proximity and dynamics of the larval 
habitats. Risk 1s further weighted by the dispersion 
and survival rates of adult mosquitoes, the 
availability of  reservoir and the human 
vulnerability (which depends on availability of 
protection and preventive devices as well as 
acquired individual immunity). 


It 1s known that environment, climate (mainly 
rainfall amount and distribution, and temperature) 
and human activities play an important role 
in  determining the vector distribution and 
vector-borne diseases epidemiology. A better 
understanding of the contribution of those 
ecological factors on the spatio-temporal dynamics 
of mosquitoes associated with malaria and dengue 
may lead to improvement of controlling measures. 
Indeed limited resources imply to focus the actions 
in places and time where they are the most useful. It 
has been argued that remote sensing technology 
has become a pre-requisite tool to assess the 
malaria burden, to model its spatio-temporal 
distribution, and plan effectively malaria control. 
The latter are key elements to implement within an 
operational system and thus to facilitate overall 
real-time monitoring of human health [Machault et 
al. 011)1. 

Along the history of tropical medicine, maps had 
held different roles that depended on the era and the 
technical possibilities. Their objectives may be: 
communication, illustration, medical research or 
tool for public health policies. From the 90”s, 
disease mapping strongly evolved with the 
introduction of the Geographic Information 
Systems (GIS) [Birchenall (2010)] that allow the 
collection, storage, Integration, analysis, and 
display of spatially referenced data. Risk maps 
have been defined in [Kitron (2000)] as “outcomes 
of models of disease transmission based on spatial 
and temporal data”, incorporating “to varying 
degrees, epidemiological, entomological, climatic 
and environmental information”, and they have 
been applied to numerous diseases for mapping a 
current situation or even anticipating outbreaks 
with Early Warning Systems [Bergquist (2001)]. 
Numerous reviews have been published in the last 
years, attempting to describe spatial and/or 
temporal mapping and modeling of vector diseases 
using GÍS and Remote Sensing (RS) information 
[Beck et al. (2000); Bergquist (2001); Hay et al. 
(2000); Machault et al. (2011)] with some inputs on 
the operational orientation of disease mapping. 

In the present study, 1t 1s thought that the conceptual 
approach (CA) of Tele-epidemiology [Marechal et 
al. (2008)] could be applied to spatio-temporal 
mapping of entomological malaria risk in urban 
settings and to dengue risk. This CA has been 
developed and patented by the French Spatial 
Agency (Centre National d'Etudes Spatiales - 
CNES) and its partners, and consists in 
monitoring and studying human and animal 
disease dynamics which are closely related to 
climate and environment variability: 1) choice of 
appropriate satellite data and dynamical models are 
evaluated, 11) all results are assessed and double- 
checked with extensive in-situ measurements 


allowing for identification of key biological 
processes. The objective concerning malaria 
was to develop a  robust operational 
methodology to produce entomological malaria 
predictive risk maps. The whole methodology 
relied on predicted vector risk areas (obtained from 
Human Biting Rate, or HBR) which were centred 
in the vicinity of predicted Anopheles gambiae s.l. 
larval habitats. The objective concerning dengue 
was to start to investigate how the CA could 
provide information on the distribution and type of 
Aedes aegypti breeding sites in a dengue epidemic 
area. 


2. METHODOLOGY 


2.1 Malaria 

The city of Dakar is the capital of Senegal and is 
located in the Cap-Vert Peninsula. The sahelian 
climate 1s modified by the proximity of the 
Atlantic Ocean, and the summer rainy season lasts 
from June to November, with average temperatures 
between 24*C and 30%C and average rainfalls of 
about 400 mm, while the cool and drier season 
lasts from December to May with average 
temperatures between 19*C and 25*C. 
Entomological field studies have been undertaken 
during September-October of 2007 and during July 
2008 to June 2010 and are detailed in [Gadiaga 
et al. (2011); Machault et al. (2009)]. Briefly, the 
monitoring study was conducted in Dakar, Pikine, 
Thiaroye and Guediawaye, in a total of 43 zones 
each covering a 200 x 200 m approximate area. 
Adult mosquitoes sampling was  carried-out 
approximately in the centre of the 200 x 200 m 
zones, once every two weeks in each studied zone 
through human landing catch, for a grand total of 
3,096 person-nights of capture, allowing to 
calculate the Human Biting Rate (HBR, number of 
Anopheles bites per person per night). Exhaustive 
search of all type of water bodies was undertaken 
every 10 days and for each zone. The latter were 
mapped using a Global Positioning System (GPS) 
device (5 m minimum precision) while physical, 
chemical and biological parameters (size, shade, 
floating and surrounding  vegetation, water 
temperature, persistence along weeks) were 
recorded, including Anopheles larvae density. 
SPOT-5 (Satellite Pour 1"'Observation de la Terre) 
images were programmed in order to be 
contemporary to the field work and were acquired 
during the three summer rainy seasons for the 
following dates: September 26, 2007, September 
24, 2008, September 28, 2009, as well as during 
one dry season (May 11, 2009). A Digital 
Elevation Model (DEM) at 90 m spatial resolution 
was available from the Shuttle Radar Topography 


Mission (SRTM versión 4.1). Decadal Moderate 
Resolution Imaging Spectroradiometer (MODIS- 
Terra) images were available at 1 km spatial 
resolution for the full duration of the field work. 
Ground measurements of daily rainfall event were 
avallable from the Senegalese weather agency for 
years 2007, 2008 and 2009. 

A GIS was built in ArcGIS 9.2 (Environmental 
Research Systems Institute, Redlands, CA, USA) 
and statistical analyses were performed using Stata 
11 (Stata Corporation, College Station, Texas). All 
ground information (water bodies and larval 
habitats, HBR, rainfall) was included in the GIS 
whilst processed images (indicators, LULC, LST, 
altitude) were added as georeferenced layers. 
Regressions between environmental and 
meteorological variables and presence of water or 
not (Step 1, logistic regression), presence of 
Anopheles gambiae s.l. larvae or not (Step 2, 
logistic regression) and An. gambiae s.l. HBR 
(Step 3, negative binomial regression), were 
fitted using all environmental and meteorological 
indicators as explanatory variables. 


2.2 Dengue 

The studied area for dengue work has been chosen 
in Tartane, a location of the Caravelle peninsula, in 
the North East of Martinique. The centre of the 
village 1s located by the seacoast but some houses 
have been built in the hills. Dengue epidemics 
that occurred in 2001 and 2010 have started in 
this area of the island. Dry season occurs from 
February to March while rainy season lasts from 
July to October. The study period was chosen from 
June 2009 to August 2011, as an extended period of 
the 2010 epidemics. 

Entomological field data consisted in the routine 
surveillance data collection undertaken by the 
“Centre de Démoustication/Lutte Anti-Vectorielle” 
of Fort de France. Briefly, field workers 
repeatedly visited houses of the studied area and 
recorded the  presence/absence of potential 
breeding sites that were filled with water, the 
number of water bodies that were positives with 
Aedes larvae and the type of water bodies. 
Regarding the known ecology of Aedes aegypti 
which is the dengue vector in Martinique and that 
1s mainly defined by small breeding sites and low 
dispersal, the experimental study unit has been 
chosen as: the house and its close environment 
(garden, yard, neighbors). Indeed, even the 
highest spatial resolution image could not allow to 
directly map the Aedes larval habitats. In 
consequence, remotely-sensed data has been 
acquired and treated in order to characterize this 
study unit and further investigate the 
environmental and meteorological determinants of 
the presence of breeding sites. The exact size of this 


study unit will further be refined using the results 
of the statistical analysis and several scales are 
expected to play a role in the presence of the 
vector (e.g. fine scale = type of house, trees in the 
garden, and large scale = density of houses in the 
neighborhood, presence of nearby vegetation). 

As the size of the experimental study unit was 
small, a very high resolution image was chosen. A 
Goe-Eye-1 image at 0.5 m spatial resolution 
has been acquired from archive images for 
13/03/2011. It has been orthorectified using the 50 
m Digital Elevation Model of IGN (Institut 
National de  l'Information  Géographique et 
Forestiere).  Pixel-based and object based 
classification were  undertaken, as well as 
calculation of several environmental (vegetation, 
humidity, soil) indices. 

Meteorological data were available from two 
ground stations located nearby the studies zone. 
They provided daily precipitation and temperature 
data. In addition, radar data was available at 1km 
spatial resolution and hourly rainfall information. 
The latter 1s of particular interest as 1t allows taking 
into account the spatial heterogeneity of rainfall. 


3. RESULTS AND DISCUSSION 


3.1 Malaria 

First step was for detection of water bodies in 
Dakar. Each observation included in the logistic 
regression  (n=48,858) corresponded to the 
presence/absence of water, based on the maximum 
water area recorded on the ground in the 45 
studied zones during the 3 studied years, at 10 
m resolution. In the global univariate and 
multivariate analysis, the SPOT-5 Modified 
Normalized Difference Water Index (NDWD of 
the rainy season and the SPOT-5 NDVI of the 
dry season were positively associated with the 
presence of water while the SPOT-5 built-up area 
and the altitude from DEM were negatively 
associated. The inversion and extrapolation of the 
model allowed drawing one map for the water 
bodies for each year, at 10 m spatial resolution. 

Step 2  highlighted the environmental and 
meteorological determinants of the presence of An. 
gambiae s.l larvae or not, recorded during the 
field work. In the global  univariate and 
multivariate analysis, the SPOT-5 NDWI Mc 
Feeters and the SPOT-3 Soil Brightness Index (BI) 
of the dry season, as well as the current night 
MODIS LST and the total ground rainfall in the 
preceding 30 days were positively associated with 
the presence of Anopheles larvae in the water 
bodies. The NDWI Mc Feeters and the Soil BL, as 
both environmental explanatory variables for the 
presence of larvae, were computed, using GIS, in 


and around each water body predicted during the 
first step. For every single day of the full duration 
of the study, including daily LST and rainfalls, the 
inversion of the multivariate global model allowed 
drawing daily maps of the probabilities of presence 
of An.gambiae s.l. larvae in water bodies. Figure 1 
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provides zoom images of the predicted probability 
of presence of An. gambiae s.l. in three districts of 
Dakar, calculated at the level of the water bodies 
predicted at step 1. 
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Figure 1. Predicted probability of presence of An. gambiae s.l. larvae (step 2) for August 
22, 2009 in the water bodies predicted at step 1 


Step 3 allowed defining the relationship between 
the predicted larval productivity surrogate and the 
field measured HBR, in order to generate 
An.gambiae s.l. HBR maps In the global 
univariate and  multivariate analysis, the 
Anopheles larval productivity surrogate (sum of 
the probability of presence of larvae calculated at 
step 2 from 30 to 1 days before the actual adult 
catching night, since the daily vector survival 
rate 1s estimated at 82% in the Dakar peninsula 


[Wercruysse (1985)], summed-up in 200 m buffers 
and rings going from 300 m to 1,000 m of the 
catching points, the weight decreasing with 
distance) and the ground rainfall in the preceding 
7 days were positively associated with the HRB 
whereas the SPOT-3 built-up and asphalt surface 
were a protective factor. The inversion of the 
global model allowed the daily prediction of the 
An.gambiae s.l. HBR at any location in Dakar with 
10 m spatial resolution (Figure 2). 
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Figure 2. Predicted An.gambiae s.l. HBR (step 3, number of bites per person per night) 
for September 22, 2009 


Modelling malaria risk requires as a first step, to 
model the entomological risk. In the present 
study, a step-by-step mapping approach has been 
followed in order to closely relate modelling and 
the biological and physical mechanisms of the 
presence of water bodies, the larval development 
and of the adult survival and dispersion. The 
temporal resolution of larvae and adult Anopheles 
mapping was daily, allowing weekly, monthly 
and/or yearly maps production that could help in 
mosquitoes control within Dakar. 


3.2 Dengue 

Photo-interpretation of the results of pixel-based 
and object-based classifications and calculation 
of kappa coefficients allowed identifying the 
strengths and weaknesses of both techniques, in 
the context of the present study. Indeed, some 
elements of the land cover had to be particularly 
well classified in the aim of characterizing the 
experimental study units (e.g. building, gardens, 
and trees). The buildings were better mapped 
using the object-based method while the accuracy 
of the classification of the vegetation, trees, and 
bare soils was higher using the  pixel-based 
method. The creation of a composite 
classification image allowed improving the global 
accuracy of the land cover classification. This 


layer has been added in a GIS containing the 
entomological data in order to further apply the 
Tele- epidemiology CA. 


4. CONCLUSION 


The Tele-epidemiology CA has been 
successfully applied for mapping entomological 
malaria risk in  sub-Saharan urban setting. 
Remotely-sensed environmental and meteorological 
data, associated with a large amount of ground 
entomological data collected specifically, allowed 
setting-up a robust operational methodology to 
draw different levels of malaria entomological 
dynamic risk maps that could be of interest for 
planning and targeting malaria 

control in urban settings. The first mapping 
level predicted the potential An. gambiae s.l. 
larval habitats in Dakar (depending upon the 
predicted presence of water, as well as 
environmental and meteorological conditions) that 
could guide local larval management. The second 
level was a map of the An. gambiae s.l. HBR 
(driven by surfaces and potential productivity of 
larval habitats as well as urbanization levels and 
rainfall amount) which could guide adult 
mosquitoes control activities and could help 


concentrating and  reinmforcing awareness and 
information messages to local inhabitants. The 
results of the present study could be seen as 
providing the basic elements for real-time 
monitoring of human health. 

Basic information has been prepared in order to 
apply the Tele-epidemiology CA to the mapping 
of dengue vectors. The expected results should 
allow to better understand the impact of 
environmental and meteorological determinants on 
the spatial and temporal distribution of Aedes. 
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RESUMEN 


Hoy en día más del 75% de la población vive 
en las ciudades, razón por la cual las áreas 
urbanas son sistemas en constante cambio que 
requieren actualizaciones periódicas sobre los 
nuevos desarrollos y el cambio de uso de suelo 
entre Otras problemáticas de los recursos 
naturales urbanos, que ayuden a generar las 
bases para las políticas públicas relacionadas 
con el desarrollo urbano sostenible. El objetivo 
principal de este trabajo fue desarrollar un 
procedimiento metodológico basado en el uso 
de sistemas de información geográfica y 
percepción remota que ayudarán a demostrar la 
utilidad de dos sensores de alta resolución 
(LIDAR y World View 2) en la cuantificación 
de los recursos naturales y urbanos en el 
Noreste de México. Los resultados forman 
parte preliminar de un esfuerzo que se 
desarrolla en el Laboratorio de Sistemas de 
Información Georreferrenciada del ITESM, en 
el desarrollo de nuevos métodos que 
establezcan mejorasen el proceso de 
actualización de información urbana y 
periurbana. Se eligieron 2 áreas de estudio de 
6.25km” en 2 Municipios pertenecientes al 
Área Metropolitana de Monterrey, San Pedro 
Garza García y Monterrey. La información 
LIDAR, consistente en una nube de puntos con 
resolución horizontal de 0.70 m y precisión 
vertical de 0.10 m de Diciembre de 2010, 
primero, se clasificó la información LIDAR y 
con ella se generó información en 3D de la 
superficie terrestre que fue clasificada. A la par 
se llevó a cabo una clasificación no supervisada 
de una imagen de satélite World View 2 
coincidente con las fechas de los datos LIDAR. 


Ambos resultados fueron validados y 
comparados a través de una matriz de error en 
la que el porcentaje de exactitud alcanzando una 
exactitud del 89%. 

Palabras clave: Datos de alta resolución, 
LIDAR, Evaluación de recursos naturales y 
urbanos 


ABSTRACT 


Nowadays over 73% of the population lives in 
the cities, reason why urban areas are 
constantly changing systems requiring regular 
updates on new developments and changes in 
land use and other issues of urban natural 
resources that help create the basis for public 
policies related to sustainable urban 
development. The main objective of this paper 
was to develop a methodological approach 
based on the use of GIS and remote sensing that 
help demonstrate the utility of two high 
resolution sensors (LIDAR and World View 2) 
for quantifying natural and urban resources in 
the Northeastern Mexico. The preliminary 
results are part of an effort that develops the 
Georreferenced Information Systems 
Laboratory ITESM, on developing new 
methods to establish improvements in the 
process of updating information urban and per 
urban. We chose two study áreas 6.25 km” in 
2 Municipalities belonging to the Metropolitan 
Area of Monterrey, San Pedro Garza Garcia 
and Monterrey. LIDAR information, consisting 
of a cloud of points with horizontal resolution 
of 0.70 m and 0.10 m vertical accuracy of 
December 2010, the first, was classified 
information with her LIDAR and 3D 


information is generated from the  earth's 
surface that was classified. The pair 
conducted an unsupervised classification of a 
satellite image World View 2 dates coinciding 
with LIDAR data. Both results were validated 
and compared through an error matrix wherein 
the percentage of accuracy reaching an accuracy 
of 89%. 

Keywords: High Resolution Data; LIDAR; 
Urban and Natural Resources Survey 


1. INTRODUCCIÓN 


Se requiere de un estudio completo sobre las 
problemáticas urbanas para afrontar los retos 
futuros tanto ambientales, económicos y 
sociales. A nivel mundial, América Latina 
reporta cerca del 3% en crecimiento urbano 
(Benítez-Senteno, 1992). El crecimiento de una 
ciudad se traduce a cambios en la superficie 
del terreno, por lo que seguramente cada mes 
o incluso cada semana existirán cambios que 
reportar por modificaciones de los elementos 
en el mapa de una ciudad; elementos tales 
como edificios, puentes, tipo de suelo, área 
pavimentada, árboles, etc.. Sin embargo, dicha 
información no es evaluada ni monitoreada 
periódicamente, complicando la actualización 
de los datos. 

El ambiente urbano sin duda es uno de los 
más cambiantes, lo cual aumenta la 
complejidad de clasificar y mapear las ciudades 
con métodos tradicionales (Williamson 
Fourie, 1998). Oficialmente el uso de suelo que 
se reporta para en las ciudades para edificios 
por ejemplo, consiste en una combinación que 
varía de áreas habitacionales (residenciales), 
comerciales oO las exclusivas del sector 
industrial. La clasificación de las áreas que 
corresponden a cada tipo de uso de suelo 
depende de datos generados a escalas 
1:50,000 ó 1:20,000 que no son tan actuales y 
que en algunas ocasiones no se encuentran en 
formatos digitales, para las ciudades algunas 
capas alcanzan 1:5,000 pero son escasas. 
Desafortunadamente, algunos de estos datos 
además presentan deficiencias en la 
información oO tienen campos o formatos 
distintos que no permiten su comparación, 
imposibilitando en algunos casos la continuidad 
del monitoreo, un problema común cuando se 
requiere hacer análisis históricos. 

La ventaja de tener resultados de análisis 
geoespaciales de la superficie, es su utilidad 
para la planeación estratégica de la ciudad 
(Chen et. al. 2008). En la mayoría de las 
ciudades latinoamericanas existen problemas 


importantes de planeación urbana, desorden de 
crecimiento, falta de regulación en los 
asentamientos, sus características de 
construcción y la rapidez con la que se 
modifica la cobertura del (CEPAL, 2001; 
PNUMA, 2003; UNCHS, 1996). Las 
consecuencias de la falta de información y de 
instrumentos legales que regulen el desarrollo 
urbano inadecuado son la pérdida de 
suelo, disminución de vegetación, 
disminución de la infiltración natural, exceso 
en los materiales de arrastre, cambios en el 
clima local, problemas de inundación, aumento 
de contaminación atmosférica (WHO, 2005), 
entre otros; todo debido a la ausencia de orden 
en el desarrollo de las metrópolis. 

Además, el análisis espacial permite evaluar 
periódicamente las bases de datos sobre el 
cambio de uso de suelo de las áreas urbanas, al 
poner en números y en planos la 
geoinformación tal como reacomodo catastral, 
áreas verdes urbanas, asentamientos irregulares 
u obstrucción a cauces, entre otras 
problemáticas de los recursos urbanos que 
permitirán guiar a los tomadores de decisiones 
con las bases necesarias para reajustar políticas 
urbanas. 

Es por lo anterior que la metodología para el 
análisis de la superficie (uso de suelo) en las 
ciudades es trascendental y requiere de un 
sistema de inteligencia que integre las variables 
de infraestructura urbana y recursos naturales en 
una base  georreferenciada que permita 
monitorear dichos cambios de una manera 
eficiente y en un costo razonable. Los sistemas 
inteligentes ayudan a conocer con datos más 
precisos el comportamiento de la estructura de 
ciudad y su proceso paralelo de transformación 
del hábitat (Kaiser et. al., 1995; World Bank, 
1995), pues permiten el análisis cuantitativo de 
los elementos de la superficie. 

Existen métodos para el mapeo de las ciudades 
tanto con las técnicas convencionales, como con 
percepción remota como lo son las estaciones 
totales, pictometría o fotografía oblicua, 
interpretación y clasificación de 1mágenes 
satelitales o más recientemente el análisis de 
datos LIDAR (Ackermann, 1999; Flood, 2001). 
LIDAR es una herramienta tecnológica para la 
creación de información 3D que está 
revolucionando desde hace una década a 
distintas disciplinas, demostrando utilidad y 
eficacia para resolver los problemas de 
precisión de datos y para el inventario de 
recursos múltiples de la ciudad (Reutebuch 
et.al., 2001; Satale Kulkarni, 2006). El 
sistema LIDAR funciona en muchas formas de 
manera más eficiente para obtener cartografía 


que la fotogrametría convencional (Baltsavias, 
1999) y los datos (nubes de puntos en formatos 
Xyz) se pueden combinar con la información de 
otro tipo de sensores. Las fusiones de la datos 
LIDAR con otros datos multiespectrales pueden 
ser aplicables a una amplia gama disciplinas 
como ingeniería civil, planificación urbana, 
manejo forestal, hidrología e hidráulica, etc 
(Filin et.al., 2005; Lefsky et. al., 2002). 

El objetivo principal de este estudio fue utilizar 
ambas tecnologías (datos LIDAR e imágenes 
multiespectrales de alta resolución) y comparar 
los resultados de la clasificación para mejorar el 
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procedimiento de obtención de información de 
la superficie de las ciudades en México. 


2. AREA DE ESTUDIO 


El área de estudio se encuentra en la ciudad de 
Monterrey, Nuevo León, México. En este 
estudio se seleccionaron 2 de los 10 sitios que 
serán estudiados utilizando datos de alta 
resolución, LIDAR y ópticos. En la Figura 1 se 
aprecian los dos sitios que numéricamente 
corresponden a los sitios 1 y 3 que tienen un 
área de 2.5 km por lado. 





Figura 1. Localización del área de estudio, los sitios de muestreo y los sitios de los que se presentan 
resultados preliminares 


3. METODOLOGÍA 


El análisis de la superficie fue realizada con dos 


metodologías distintas de percepción remota: 


1) Clasificación de datos LIDAR y 2) 
Clasificación de imágenes ópticas (World View 
2). La metodología LIDAR fue la siguiente: Para 
el filtrado de los datos crudos se utilizaron 


algoritmos de filtrado que tiene el 
programa  MARSs de Merrick. Las 
breakinglines permitieron separar el suelo 
(Digital Terrain Model, DTM) del resto de los 
elementos de la superficie (Digital Surface Mode, 
DSM). El resto de los elementos del DSM se 
clasificó diseñando filtros compuestos con los 
siguientes algoritmos: Building, Elevation, 
Heightfromsurface,  Intensity,  PlaneAnalysis, 
Reclassify, Return y Cross Section. Las clases 
establecidas fueron: Edificios, Calles y 
avenidas, Vegetación, Suelo, Agua, y Lecho del 
río etc. El diseño de los filtros y sus resultados 
específicos pueden ser consultados en la 
publicación de los autores Yépez Rincón é 
Lozano García, 2012, in pub. La clasificación 
de las imágenes World View 2, siguió el 
procedimiento estándar de PCI Geomática para la 
clasificación de imágenes ópticas. Realizándose 
un manejo previo de la información en la que se 
hizo la corrección geométrica y se extrajo el 
polígono con los sitios de estudio específicos. 


Para la clasificación se utilizaron 60 clases 
para hacer los clusters y con esquema de color 
en verdaderos y 24 iteraciones, y convergence 
threshold 0.95. La imagen que se manejó 
contenía únicamente cuatro de las ocho bandas 
disponibles: Roja, Verde, Azul e  Infraroja 
cercana. Las clases obtenidas fueron a) Edificios, 
b) Vegetación, 3) Calles, 4) Suelos, 5) 
Agua, 6) Sombras y 7) Otros como se 
desglosan en la Figura 2. La clasificación se 
validó al seleccionar 10 puntos por cada tipo de 
clase, es decir una estratificación de puntos al 
azar. Se realizó la misma metodología para la 
clasificación de World View 2 y para LIDAR. 
Dado a que el objetivo principal de este 
documento fue utilizar ambas tecnologías y 
comparar los resultados de la clasificación se 
utilizaron una serie de puntos estratificados por 
clase al azar. Cada clase obtuvo 10 puntos que 
fueron muestreados al azar para verificar su 
certeza. 


Clasificación y análisis de la superficie 


LIDAR 


WorldView 2 


Figura 2. Metodología que se siguió para clasificar la información 


La verificación se realizó en esta etapa sobre La información se comparó en una matriz de 


la imagen, sin embargo para el resultado final confusión donde se obtuvieron los resultados 
del estudio se pretende verificar en campo al de las ópticas en negro y de LIDAR en rojo 
menos el 50% de los puntos señalados. (Figura 3). 
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Figura 3. Las clasificaciones por sensor y sitio 


4. RESULTADOS Y CONCLUSIONES 


Ambas clasificaciones fueron comparadas 
por clase, pero para el cálculo de errores de 
omisión y comisión se generalizaron las clases 
en 1)Edificios, 2) Vegetación urbana, 3)Calles y 
avenidas, 4)Suelo y 5) Agua y 6) Otros. Se 
obtuvieron diferencias importantes entre ambos 
sensores. Desde la vista aérea se aprecia en 
la Figura 3 una clara confusión del filtrado 
en LIDAR con las clases de vegetación y 
edificios. La razón de esta confusión se remite 
al número de retornos que se utilizaron y quizás 
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a valores de intensidad. Sin embargo, la 
clasificación de la imagen se pueden apreciar 
las diferencias entre 22 clases, lográndose 
apreciar claramente entre clases con edificios 
residenciales e industriales. 

En la Figura 4 se observa con puntos rojos los 
errores del sensor LIDAR y con puntos negros 
las imágenes ópticas. Para los resultados de 
LIDAR puede concluirse que la vegetación 
tiende a confundirse con los edificios, así como 
los pavimentos con suelos. 
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Figura 4. Matriz de error para LIDAR (puntos rojos) y World View 2 (puntos negros) 


5. RECOMENDACIONES 


Los resultados preliminares obtenidos de 
acuerdo al análisis de la matriz de errores 
de omisión y comisión indican que será 
necesario el perfeccionamiento de los filtros 
utilizados en la clasificación de los datos 
LIDAR, para lograr mejores resultados en el 
filtrado automático de la información. En 
términos generales se recomienda en primer 
término 1) Ajustar la clasificación de la imagen 
World View para alcanzar un nivel de 
exactitud sobre el 95%, mediante un proceso de 
clasificación mixta (método supervisado y no 
supervisado, 2) Recalibrar los filtros que fueron 
diseñados para la clasificación de datos LIDAR 
en base a las clases propuestas para los datos 
Ópticos y 3)  Fusionar la información 
resultante de la clasificación para elevar el 
nivel de exactitud al 98%. 
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RESUMEN 


El Parque Nacional los Glaciares se encuentra 
ubicado en la región conocida como Andes 
Australes del territorio argentino, al sudoeste de 
la provincia de Santa Cruz, en el límite con 
Chile, el mismo fue declarado "Patrimonio 
Mundial" por la UNESCO en el año 1981. Gran 
parte de este Parque estuvo cubierto por 
glaciares y los ríos de hielo dieron forma al 
paisaje con montañas de laderas abruptas y 
amplios valles. El área de estudio corresponde 
al Glaciar Perito Moreno, ubicado sobre el 
Brazo Sur del Lago Argentino, al Glaciar 
Upsala, ubicado sobre el Brazo Norte del mismo 
lago y a las áreas aledañas a los mismos. El 
uso de geotecnologías basado en la utilización 
de información satelitaria representa una 
alternativa de Importancia práctica como 
herramienta de monitoreo y evaluación de 
recursos naturales. El objetivo de este estudio 
fue realizar un mapeo temático de los recursos 
naturales presentes en el área, para identificar 
y caracterizar los tipos de ambientes y las 
unidades  fisonómicas forestales nativas, 
generando una cartografía temática de las 
unidades identificadas. Para generar dicha 
cartografía se utilizó información 
correspondiente a imágenes radar (RADARSAT 
2,  polarimétricas), en combinación con 
imágenes ópticas Landsat 5 TM; se utilizaron 
técnicas de procesamiento digital, 
clasificaciones  multiespectrales y análisis 
polarimétricos, integrando los datos de radar y 
Ópticos, aprovechando el aporte que ofrecen 
ambas tecnologías. La cartografía temática 
generada será integrada en un SIG territorial. 
Palabras claves: SAR / cartografía / bosques / 
geotecnologías / SIG 


ABSTRACT 


Los Glaciares Argentina National Park is 
located in Andes Patagonicos and it belongs to 
Santa Cruz province and it limits with Chile 
country. It was declared Natural Patrimony of 
Mankind by the United Nations Educational 
Scientific and Cultural Organization (UNESCO) 
in 198l. An important part of this Park is 
covered by glaciers and the ice rivers give form 
to the landscape with abrupt hillsides mountains 
and wide valleys. Glacier Perito Moreno 
belongs to the study area and it is located over 
Brazo Sur that belongs to the Lago Argentino. 
There 1s another study area Glaciar Upsala 
located over Brazo Norte of the same lake. 
Besides, the environment of both glaciers 1s 
taken into account to carry out a complete study 
of this important area. Satellite data has probed 
to be very important for monitoring and 
evaluating natural resources. The objective of 
this paper 1s to carry out a thematic map of the 
natural resources that belong to this area. In this 
way the different environment can be identified, 
characterized and the physiognomic natural 
forest unities described so a  thematic 
cartography of the different identified unities 
are obtained. Polarimetric RADARSAT 2 and 
optical Landsat 5 TTM images were used to 
generate this cartography. Digital processing, 
supervised classification, polarimetric analysis 
techniques are applied to integrate radar and 
optical data, taking into account both 
technologies at the same time. All the thematic 
cartography will be added to a GIS. 

Keywords: SAR / cartography / forestry / 
geotechnologies / GIS 


INTRODUCCIÓN 


La necesidad de obtener información sobre 
recursos naturales de carácter forestal y de 
ambientes asociados, es fundamental para 
establecer políticas de acción que optimicen las 
funciones y actividades relativas a la 
conservación y el aprovechamiento económico 
de tales recursos, así como a su control y 
fomento, logrando de esta manera un manejo 
sostenible de los mismos. 

La teledetección espacial se ha mostrado como 
una herramienta de importancia práctica para la 
evaluación y el monitoreo de los recursos 
forestales, acrecentada por el desarrollo 
tecnológico de los últimos años que ha 
posibilitado una amplia disponibilidad de 
información de distintos tipos de sensores, 
procesos y algoritmos que permiten obtener una 
gran variedad de productos utilizados en 
proyectos de gran complejidad. 

Este trabajo está orientado a la utilización de 
imágenes ópticas y de radar para obtener 
información temática de cubiertas naturales; en 
el caso de las aplicaciones forestales el uso de 
ambos datos, ópticos y de radar se hace 
necesario ya que permite superar limitaciones 
que cada uno de los tipos de datos tienen por 
separado; aunque existen métodos de obtención 
de índices y clasificación de vegetación 
mediante datos ópticos, estos no pueden ser 
utilizados en regiones cubiertas frecuentemente 
por nubes como lo es la región andino- 
patagónica. La señal radar retrodispersada por 
una superficie vegetal depende de dos tipos de 
factores: (1) factores geométricos relativos a la 
estructura del suelo y de la cubierta vegetal y, 
(1) factores  dieléctricos principalmente 
controlados por la humedad del suelo y del 
dosel (Dobson, 1995). En términos generales, 
las superficies más secas producirán una señal 
más baja. La influencia de las propiedades 
estructurales sobre la señal devuelta al sensor 
por una cubierta vegetal es mucho más 
compleja puesto que se ve afectada por la 
frecuencia, la polarización y el ángulo incidente 
de la señal enviada por el sensor. 

Los parámetros estructurales de mayor 
importancia son la arquitectura del dosel 
(tamaño, forma y orientación de elementos 
retrodispersantes: hojas, ramas y troncos) y la 
rugosidad del suelo. Por otro lado la topografía 
constituye otro factor que causa dificultades 
para identificar los distintos estratos 
altitudinales de vegetación y que debe se 
rconsiderada en el momento del análisis de los 
resultados. En este sentido, este trabajo 
demuestra la complementariedad de ambos 


tipos de datos en estudios de cubiertas 
forestales. 

En este contexto, se plantea como objetivo 
general del presente trabajo generar una 
cartografía temática de las unidades de 
vegetación existentes en el área de estudio que 
puedan ser identificadas y  categorizadas 
mediante el análisis digital y reconocimiento a 
campo de los diferentes tipos forestales: 
(Nothofagus  pumilio,  lenga;  Nothofagus 
Antarctica, ñire y Nothofagus  betuloides, 
guindo); la vegetación de estepa (arbustos 
pequeños, espinosos y pastizales) y una zona de 
transición (especies típicas del bosque y de la 
estepa patagónica). Asimismo es de interés el 
mapeo de otras cubiertas presentes en el área 
tales como suelo desnudo (roca), hielo de 
glaciar (zona de ablación y acumulación) y 
nieve, con el objetivo de evaluar los 
mecanismos de interacción de estas cubiertas 
con la señal del radar polarimétrico, utilizando 
metodologías de procesamiento y análisis de 
imágenes. El desarrollo del mismo permitirá 
generar información para estudios científicos y 
técnicos, así como para la planificación e 
implementación de estrategias de gestión 
sustentable del Parque Nacional, coherentes con 
sus características ambientales. 


MATERIALES 


Área de estudio 

El área de estudio, Parque Nacional los 
Glaciares, Provincia de Santa Cruz, Argentina 
se localiza entre los 49% 43” y 50” 53” de latitud 
Sur y los 71% 52” y 73" 43” de longitud O 
(Figura 1). 

En ella se pueden distinguir fisiográficamente 
tres fajas que corren en forma paralela al eje de 
la cordillera de sur a norte: 1- Cordillera 
Principal o región andina, 2- Región Pre-andina, 
3- Región Sub-andina. 

La primera corresponde a una extensa masa 
montañosa con alturas superiores a los 2500 m, 
su mayor ancho se estima entre 50 y 60 km. En 
la parte central está cubierta por campos de 
hielos denominado Hielo Continental 
Patagónico, de los cuales descienden a dos 
grandes lagos, el Lago Argentino y el Lago 
Viedma, y lenguas glaciares, entre las cuales 
podemos identificar 13 grandes glaciares. 

La segunda región está conformada por un gran 
valle que desciende desde el cerro Norte hasta la 
estancia Cristina, el brazo Norte, el brazo 
Spegazzini, el Canal de los Témpanos y el brazo 
Sur (todos ellos brazos del Lago Argentino). 


Las cadenas montañosas bajan desde alturas 
superiores de los 2000 m y se encuentran 
desprovistas de hielo. 

La tercera región sub-andina se compone de 
mesetas que se levantan abruptamente sobre el 
nivel del Lago Argentino y van disminuyendo 
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hacia el este confundiéndose gradualmente con 
la altiplanicie patagónica (Feruglio, 1944). 

En esta etapa del proyecto se seleccionaron 
como áreas pilotos: el Glaciar Perito Moreno, el 
Glaciar Upsala y zonas aledañas a ambos 
glaciares. 





Figura 1. Localización del área de estudio 


Imagen óptica 


La imagen óptica seleccionada fue la del 
Landsat 5 TM, Path/Row 231/095, de fecha 20 
de febrero de 2011, con ancho de barrido de 
185 km y resolución espacial de 30 m, del 
banco de datos del Centro de Sensores Remotos, 


provista por CONAE (Figura 2), la misma fue 
escogida considerando características de baja 
nubosidad, buena calidad radiométrica y con un 
intervalo de tiempo entre las imágenes relevante 
al problema en estudio. 





Figura 2. Imagen Landsat 5 TM FCC (R/Irm/Irc — A/V/R) 


Imágenes radar 


Las imágenes empleadas en el presente estudio 
fueron adquiridas por el sensor SARa bordo 
del satélite RADARSAT-2. Se contó con dos 
imágenes RADARSAT-2, SAR, ID 36726 e 
ID 36757, ambas del 27 de abril de 2009, 
adquiridas bajo el programa Science and 
Aplications Research (SOAR) de la Agencia 
Espacial Canadiense. El satélite RADARSAT-2 
fue lanzado en diciembre de 2007 y está 
gestionado por la empresa MacDonald- 
Dettwiler and Associates Ltd. (MDA) y el 
gobierno de Canadá. 

El tipo de imagen empleada en este trabajo es el 
denominado Fine Quad-Polarization con cuatro 
canales de polarización (HH, VV, HV, VH), 
ancho de barrido nominal 25 km y una 
resolución espacial de 10 m. 

El SAR  interferométrico y  polarimétrico 
(PolInSAR) utiliza, además, la información 
proveniente del hecho de que la retrodispersión 
es una onda electromagnética que tiene un 
ángulo preferente de polarización, dependiendo 
de las características de la superficie. Utilizando 
filtros de polarización tanto en la emisión como 
en la recepción de la señal, es posible estudiar la 
misma, en cuatro configuraciones angulares 
para cada punto de la superficie: HH, VV, HV y 
VH. 

El PolInSAR utiliza una información muy rica, 
pues para cada punto de la superficie dispone de 
información no sólo del módulo de la 
retrodispersión, sino también de su fase y de sus 
características angulares. Con esta información 
es posible procesar los datos basándose en 
modelos físicos que predicen la retrodispersión 
en función de una serie de parámetros de la 
superficie que pueden ser medidos en el terreno. 


Software 


- ENVI 4.7 

- NEST 4.C (ESA) 

- POLSARPRO V4.2.0 (ESA) 
- ARC-VIEW 3.2 


METODOLOGÍA 


Para caracterizar el área de estudio, se realizó 
una revisión de los trabajos existentes 
relacionados con la utilización de datos ópticos 
y de radar en ambientes de bosques y glaciares, 
su cartografía e interpretaciones. A partir de 
esos datos y del análisis de imágenes satelitales 
provenientes de diferentes sensores, Ópticos y 
radar, se estructuró una nueva serie de mapas 
temáticos correspondientes a las comunidades 
de vegetación y de ambientes presentes. 

Se evaluaron a través del análisis visual, sitios 
característicos de los diferentes ambientes 
presentes en el área de estudio utilizando como 
apoyo para este trabajo productos mejorados a 
través de procesos digitales simplificados 
(Ensanche de contraste y Realce de bordes), 
aplicados a la imagen Landsat 5 TM empleada. 
Se realizó un relevamiento de campo con el 
objetivo de identificar y validar cada una de las 
clases temáticas que fueron establecidas en el 
trabajo de interpretación visual preliminar y 
antecedentes de trabajos realizados en el área de 
estudio, el trabajo de terreno consistió en recorrer 
el área de estudio, localizando con GPS muestras 
de los distintos tipos de cobertura vegetal boscosa 
(Bosque de lenga, Bosque de ñire y Bosque 
asociado), vegetación de estepa, mallines y otras 
cubiertas asociadas como hielo de glaciar, nieve, 
roca, etc. (Figura 3). Estos sitios de entrenamiento 
sirvieron como base para definir patrones 
espectrales a ser utilizados en la clasificación 
digital supervisada de la imagen Landsat y para el 
análisis de los mecanismos de interacción en los 
datos radar. 





(d) 





(2) 


Figura 3. Fotos de relevamiento a campo. (a) Bosque de lenga, (b) Bosque de fire, (c) Bosque asociado, (d) 
Vegetación de estepa, (e) Mallines, (£) Roca y (g) Hielo de glaciar y Nieve 


La 1magen Landsat fue ajustada 
geométricamente mediante el método imagen- 
imagen, utilizando como referencia una escena 
Landsat 7 ETM +, ortogeoreferenciada, WRS-2, 
path/row 231/95, del 14 de octubre de 2001, el 
sistema de referencia corresponde a la 
proyección UTM, zona 18 sur, con Datum 
W6GS84 (Fuente USGS/ GLCH). 

Se realizó una clasificación digital 
multiespectral supervisada, considerando las 
clases temáticas: bosque de lenga, bosque de 
ñire, bosque asociado, vegetación de estepa, 
mallines, roca, hielo glaciario y nieve. Se 
verificó la buena separabilidad de clases y se 
redistribuyeron las mismas. Para la evaluación 
del comportamiento y separabilidad espectral de 
las distintas clases se aplicó el análisis de 


divergencia transformada. El algoritmo 
utilizado en la clasificación fue el de Máxima 
Verosimilitud. Las imágenes RADARSAT 2, se 
calibraron radiométricamente, por la variación 
del ángulo de incidencia. El ángulo local de 
incidencia del haz aumenta gradualmente entre 
los puntos situados en el rango cercano y el 
lejano. Asumiendo el caso más sencillo de una 
superficie plana, esto produce una atenuación en 
la respuesta para los puntos situados en el rango 
lejano. Para compensar este efecto es que las 
imágenes de radar deben ser calibradas y se le 
aplicó el filtro Frost con ventana de 3x3 para 
disminuir el ruido speckle. 

Se realizó el análisis polarimétrico, la 
descomposición de Pauli y de Sinclair en las 
imágenes RADARSAT-2. El análisis 


polarimétrico es relevante solamente si parte de 
la onda retro-dispersada es predominante. La 
evaluación de las matrices de retrodispersión y 
coherencia permite cuantificar las dos partes de 
la onda (polarizada y no polarizada) (Chan, 
1981). La descomposición de los blancos 
permite la definición de los parámetros 
polarimétricos tales como el análisis de la 
Entropía, Anisotropía y Alpha; estos parámetros 
reúnen la información polarimétrica y la 
descomponen en términos de los mecanismos de 
retrodispersión que se producen en las cubiertas, 
lo que está directamente relacionado con su 
naturaleza (Cloude, 1996). La Entropía (H) 
representa el grado de aleatoriedad o desorden 
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estadístico de la dispersión. Varía entre 0, 
reflexiones de primer orden, y 1, mezcla 
aleatoria de mecanismos de reflexión. El ángulo 
Alpha proporciona un valor relativo al 
mecanismo de retrodispersión dominante, oscila 
entre 0%, dispersión superficial y 90%, reflexión 
de esquina, con valores intermedios 45" 
ilustrando la dispersión de volumen. Por último, 
la Anisotropía (A) cuantifica la proporción entre 
los dos mecanismos de dispersión no 
dominantes. Su valor varía entre O y 1 y da idea 
de la homogeneidad de la cubierta; valores 
bajos corresponden a cubiertas homogéneas y 
viceversa (Figura 4, 5, 6 y 7). 





(c) 


Figura 4. Parámetros polarimétricos en la imagen RADARSAT-2, 
correspondiente al área del glaciar Perito Moreno. (a) Alpha, (b) Anisotropía y (c) Entropía 
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Figura 5. Parámetros polarimétricos en la imagen RADARSAT-2, 
correspondiente al área del glaciar Upsala. (a) Alpha, (b) Anisotropía y (c) Entropía 





Figura 6. Descomposición de Pauli (a) y Descomposición de Sinclair (b) en la imagen RADARSAT-2, 
correspondiente al área del glaciar Perito Moreno 
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Figura 7. Descomposición de Pauli (a) y Descomposición de Sinclair (b) en la 1magen RADARSAT-2, 
correspondiente al área del glaciar Upsala 


El teorema Pauli, utiliza las matrices de Spin 
Pauli para expresar [S] en términos de 
dispersión única (superficie plana, esfera, oO 
triedro) y mecanismos diédricos con una 
orientación de 0% y 45%, estos algoritmos 
permiten descomponer las observaciones en 
parámetros que ofrecen una interpretación física 
de la interacción ocurrida (por ejemplo 
dispersión superficial, dispersión de volumen o 
doble rebote) (Freeman, 1992; Lee, 19099; 
Ouarzeddine, 2007). Las descomposiciones se 
pueden utilizar para interpretar la naturaleza de 
las cubiertas e incluso realizar clasificaciones 
supervisadas y no supervisadas (Alberga, 2007; 
Rodrígues, 2003). 


RESULTADOS 


A partir del análisis de la cartografía vigente, la 
interpretación visual y de la clasificación digital 
multiespectral de la imagen Landsat (Figura 8 y 
Tabla 1) se generó un mapa temático con las 
categorías objeto de estudio, en este caso el 
resultado es el mapa de unidades de vegetación 
y ambientes, considerando las clases temáticas: 
bosque de lenga, bosque de ñire, bosque 
asociado, vegetación de estepa, mallines, roca, 
hielo glaciario y nieve. 

Clases temáticas identificadas: 

Bosque de lenga: Corresponde a la especie 
lenga (Vothofagus pumilio) como especie 


predominante, en dos tipos:  arbóreo y 
achaparrado. Según donde crezca la lenga y en 
función de su distribución altitudinal puede 
presentar gran porte y diámetro hasta 
transformarse, en la altura de las laderas, en un 
arbusto de pocos metros de altura. 

Bosque de ñire: Corresponde a la especie ñire 
(Nothofagus antartica), es una especie bastante 
rústica, de la familia de las fagáceas, se 
encuentran bosques puros, con árboles de hasta 
15 m de altura, en la zona de transición entre el 
bosque y la estepa y también asociado con la 
lenga en los faldeos. También puede presentarse 
en forma arbustiva. 

Bosque asociado: Corresponde a un bosque 
mixto, que crece en sitios con condiciones más 
favorables, como faldeos con mayor insolación 
y cercanos a orillas de lagos, está conformado 
principalmente por otra especie de la familia de 
las fagáceas, el guindo, (Vothofagus 
betuloides), con individuos de gran porte, 
perennifolio, alcanzando los 20 a 30 m de altura 
y asociado a la lenga,  notro, (Embothrium 
coccineum), canelo, (Drimys winteri) y sauco 
(Sambucus spp.). 

Vegetación de estepa: La clase vegetación de 
estepa se compone de especies de herbáceas 
perennes como los coirones y pastizales que 
forman matas bajas y compactas, con presencia 
de subarbustos bajos formando matas 
circulares, como el neneo (Mulinun spinosum), 





la mata mora (Senecio filaginoides), la mata 
negra (Junellia tridens) etc. 

Mallines: Son praderas húmedas que se 
encuentran en la estepa, generalmente asociadas 
con los cursos de ríos o arroyos o con los fondos 
de los valles. En ellos, la alta disponibilidad de 
agua, debida a la redistribución local, determina 
una fisonomía enteramente diferente. La 
cubierta vegetal a menudo supera el 100%, y 
las especies dominantes son los pastos 
mesofíticos (Poa  pratensis,  Deschampsia 
flexuosa, etc.), los juncos (Juncus balticus) y las 
ciperáceas (Carex spp.). 

Roca: áreas de roca pelada o suelo expuesto, se 
distinguen varios sectores rocosos: en la parte 
alta de las montañas, luego de los bosques de 
lenga achaparrada y en sectores bajos como la 
Península Herminita o alrededores de la 
Estancia Cristina donde se denotan los efectos 
de la erosión glaciaria. 

Nieve: áreas cubiertas de nieve ubicadas en la 
zona de acumulación de los glaciares. 

Hielo de glaciario: áreas de hielo y lenguas 
elaciarias. Están unidos en la zona de 
acumulación por un campo de hielo 
denominado Hielo Patagónico Sur, con una 
altitud promedio de 1600 m y un área de 13000 
km?. Esta clase identifica el hielo de los 
glaciares en la zona de ablación, incluyendo las 
morenas. 
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Figura 8. Imagen Landsat 5 TM clasificada 


Tabla 1. Resultado de la clasificación. Superficie por clase temática 


Average accuracy = 97.08 / Overall accuracy = 98.21 / KAPPA COEFFICIENT = 0.98251 





Se realizaron clasificaciones polarimétricas con 9 clases, a fin de observar distintas 
supervisadas y no supervisadas características de hielos, vegetación y rocas, en 
(entropía/alfa/anisotropía), Wishart-H/A/Alpha, ambas áreas de estudio (Figuras 9 y 10). 
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Figura 9. Clasificación Whishart. (a) no supervisada y (b) supervisada. Perito Moreno 
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CONCLUSIONES 


La clasificación supervisada realizada a la 
imagen Landsat 5 TM, arrojó buenos resultados 
utilizando las mismas clases temáticas, en forma 
visual y por las estadísticas (separabilidad de 
clases, matriz de confusión). 

En el análisis polarimétrico realizado sobre las 
imágenes RADARSAT-2 se observó mayor 
despolarización (mayores componentes de 
polarización cruzada HV, VH) en las zonas de 
bosques, según lo esperado, y también en 
algunas zonas de hielos de glaciar. 

Es imprescindible considerar y corregir los 
efectos del relieve en las imágenes SAR de 
bosques en áreas montañosas o con importantes 
gradientes topográficos. 

Si bien los mecanismos de interacción entre las 
microondas y los componentes de la vegetación 
son complejos, el tipo de información ofrecida 
por los datos SAR demuestra su importancia 
como dato único y complementario para 
estudios forestales. El desarrollo de radares en 
Banda L, polarimétricos, como el SAOCOM 
Argentino, promete aumentar la utilidad de 
estos datos para estudios de bosques y 
contribuir en aplicaciones operativas. 

Como conclusión final se demostró la 
importancia del aporte de la información radar 
aún en zonas de alto relieve con presencia de 
deformaciones geométricas, para la confección 
de mapas temáticos. 
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